00 

o 

O 
(N 






> 



S 



Analytic treatment of CRIB Quantum Memories for Light using Two-level Atoms 

J. J. Longdell 
Department of Physics, University of Otago, Dunedin, New ZealancQ 

G. Hetet and P. K. Lam 

ARC COE for Quantum- Atom Optics, Australian National University, Canberra, ACT 0200, Australia 

M. J. Sellars 
Laser Physics Centre, RSPhysSE, Australian National University, Canberra, ACT 0200, Australia 

(Dated: 1 July 2008) 

It has recently been discovered that the optical analogue of a gradient echo in an optically thick 
material could form the basis of a optical memory that is both completely efficient and noise free. 
Here we present analytical calculation showing this is the case. There is close analogy between the 
operation of the memory and an optical system with two beam splitters. We can use this analogy to 
calculate efficiencies as a function of optical depth for a number of quantum memory schemes based 
on controlled inhomogeneous broadening. In particular we show that multiple switching leads to a 
net 100% retrieval efficiency for the optical gradient echo even in the optically thin case. 



'rl'. I- INTRODUCTION 

The ability to store and recall quantum states of light as coherences in atomic media is currently being actively 
fi ' pursued. Such quantum memories would find use both in optical based quantum computation |l| and long distance 
quantum communication [^. 

Many of the current quantum memory approaches involve the use of three level atomic systems, and store quantum 
information between two quasi ground states [3, Ij, Is], la, |7| . Quantum states have been stored and recalled using this 
approach 0, S B [l3| ■ Three level schemes are particularly well suited to gaseous atomic systems where significant 
optical depths can be obtained for the optical transitions due the allowed transitions and long coherence times can be 
obtained for the ground state coherences. Rare earth ion doped solids at cryogenic temperatures have much higher 
^^ ' atom densities and allow reasonable optical thickness to be obtained from ensembles of weak oscillators. These weak 
^) , oscillators can have correspondingly long coherence times. This means that approaches involving only two level atoms 
C^ • can be considered. Sangouard et al. [ll| showed that in principle 54% efficiency can be achieved from a controlled 
reversible inhomogeneous broadening (CRIB) echo [a, |7|| with two level atoms, in the case where the broadening 
(^ ■ mechanism is 'transverse', that is not correlated with position along the optical path. The efficiency was Hmited by 
00 [ reabsorbtion, as the sample is made more optically thick, which is required to absorb the input pulse, more of the 
^^ ■ echo gets absorbed before it can make it out of the sample. Shortly afterward it was shown that the optical gradient 
Kj^ ' echo or "longitudinal" CRIB echo |12|, the only CRIB echo that has so far been reported experimentally, did not suffer 
•^H , from reabsorbtion problems and offered noise free, potentially 100% efficient storage [lJ|. An attractive property of 
p\ ' these echo based techniques is that time bandwidth product scales better with optical depth [2, [l3|- Indeed storing 
and recalling multiple pulses has been demonstrated for CRIB echo memories [l5i| but this has proved difficult for 
EIT. 

Following on from [l3| here we present analytic theory of the optical gradient echo or "longitudinal CRIB" echo. 
Analytic solutions of the Maxwell-Bloch equations are derived for a pulse of light interacting with ensemble of atoms, 
where the resonant frequency atoms varies linearly with the propagation distance. After the input pulse atomic 
coherence dephases due to the different resonant frequencies of the atomic ensemble. If the detunings of each of 
the atoms are reversed the ensemble rephases and produces an echo. In the case where the medium is optically 
thick the echo is totally efficient. Indeed, we show very close analogy between a gradient echo memory and a 
pair of beamsplitters. We use this analogy to calculate the efficiency of the longitudinal echo for optically thin 
samples. Multiple switching and a number of memories operated in series can also be easily treated using this 
beamsplitter analogy. As transverse broadening can be considered as the limit of many stacked optically think 
transverse broadening. Transverse broadening can be treated using the same approach. We reproduce previous 
results for the efficiency of transverse CRIB echoes [ll| and extend to the case of multiple switching. 
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II. THEORY 

We consider the interaction of a collection of two-level atoms with a light field where a detuning of the atoms that 
is linearly dependent on their position can be introduced. We shall assume that the area of the incoming pulses is 
much less than that of a tt pulse. This enables us to treat the atoms as harmonic oscillators. Before the detuning of 
the ions are flipped, the Maxwell-Bloch equations in the frame at the speed of light are |16| 

d 

—a{z,t) = -{-f + if]z)a{z,t)+igE{z,t) (1) 

j-E{z,t) = iNa{z,t) (2) 

oz 

Where E represents the slowly varying envelope of the optical field; a the polarization of the atoms; N the atomic 
density; g the atomic transition coupling strength; 7 the decay rate from the excited state and t]z the detuning from 
resonance. We shall assume that the process happens fast compared to the atomic decay rate and will take 7 = 0. 

After the detuning of the ions is flipped to recall the light pulse the Maxwell-Bloch equations describing the dynamics 
are the same as above except the sign of the irjz term in [l] is flipped leading to (with 7 = 0) 

— a{z,t) = —iriza{z,t) + igE{z,t) (3) 

^Eiz,t) = iNa{z,t) (4) 

As mentioned above, although the treatment here is classical, the linearity of Eqs. I|1|2|3|4I ) ensures that exactly 
the same analysis would be true for operator valued a and E. The only added noise in the output pulses will be the 
vacuum noise added to preserve commutation relations in much the way as light interacting with a beam-splitter. 
Our results will thus hold for quantum mechanical flelds also. 

The domain over which these equations will be solved is illustrated in Fig. [TJ The boundary conditions are the 
input pulse E{z — — zq, t) = f'mi^) ^^^ the initial state of the atoms. We shall assume that the atoms start in their 
ground state so that a{z, t = —00) = 0. By propagating forward in both z and t from these boundary conditions one 
could determine the flnal state of the atoms and the light exciting the sample at z = zq. Ideally the procedure would 
result in an output pulse, E{z = +zo,t) = /out(0 that has the same energy and is closely related to the input pulse, 
and that at the end of the process all the atoms are left in their ground states, a{z, t = -l-oo) = 0. It should be noted 
that we are dealing with the situation where there is no decay of the atomic states, so in order for a(z, t — +00) = 
the atoms must be left in their ground states at the end of the retrieval process, it cannot happen via decay of the 
atomic states. 

Rather than propagating forward all the way from our boundary conditions at z = — zo and t — — cxd to arrive 
at expressions for z = zq and t = -|-cxd we will take a different approach that better makes use of the symmetry of 
the situation. We flrst solve for the behavior in Region A using the boundary conditions a{z, t = —00) = and 
E{z — —ZQ,t < 0) = /in-(0- This analysis will give us expressions for E and a at the time when the detunings are 
flipped (a(z, t — 0) and E{z, t — 0)) as well an expression for the light that leaves the sample before the fleld is flipped 

(i?(z = -hzo,i<0)-/out-W)- 

We then carry out what turns out to be a very similar analysis. Working in Region B, we start with an arbitrary 
form output pulse {E{z = -\-ZQ,t > 0) = /out+(^)) ^^^ the requirement that the atoms end up in their ground states 
(a(z, t — Qo) — 0) and then we work out what is needed from the other two boundary conditions E(z = — zq, t > Q) and 
{a{z, t — Q) and E{z, t — 0)) in order for these outcomes to occur. By comparing these required boundary conditions 
with those that actually occur we show that in the case of an optically thick sample the storage is completely efficient. 
For the case of an optically thin sample we can calculate how efficient the process will be. 

The analysis of the behavior in Region A is as follows: First Eq. [T]is integrated to give 



aaiz,t) = ig dt'e-'^''^*-''^Ea{z,t') 

J — 00 

/CO 
dt' H{t - i')e""'^^*"*'^^a(2, t') (5) 

-00 

here H denotes the Heaviside step function. The above expression is in the form of a convolution, taking the Fourier 
transform gives a product. 
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Figure 1: The domain over which the Maxwell-Bloch equations are solved. Each horizontal slice represents one position in 
the sample of atoms as a function of time. The light enters the sample at z = —zo and leaves at z = +zo. The dotted line 
represents the point in time at which the detunings of the ions are flipped. In Region A the dynamics are described by Eqns. 
[Hand El in Region B the dynamics are described by Eqns. [3]and|4l 



aa(z, w) = igEa{z, uj) — ; + Trd{uj + rjz) 

\i\uj + rjz) 



(6) 



Substituting this in Eq. [2] we get 



d:,Ea{z,uj) = -g'^N[ — -^ + 7r(5(w + r/z) ] Ea{z,io) 

\i{LO + rjz) 



Integrating this equation we have 



Ea{z,ijj) = Eaiz = -zo,uj)e^p dz' -g^Nl- ■ + TrS{uj + rjz')] 

= Fin_(w) exp(-/37r(iJ(cj + r]z) - H{uj + rjzo))) 



u! + rjz 



to + rjzo 



if} 



(7) 



(8) 



where /3=^. 

It can be seen that the amplitude of each spectral component is attenuated by a factor exp(— 7r/3) after traveling 
past the position in the sample where it is resonant with the atoms. It also receives a phase shift as it travels through 
the sample. 



We make the assumption that the spectral coverage of the sample is large compared to the optical depth. That is 
for each frequency of interest lo in our input signal we have /3w <^ rjzo in which case our expression for Ea{z, cu) takes 
the form 



Ea{z,uj) = F^^_{uj) exp{-PTrH{uj + rjz)) 



uj + rjz 
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(9) 



Substituting z = +zq in the above equation one finds that the transmitted pulse is equal to the the incident pulse 
multiplied by an attenuation factor exp(— /37r). 



/out^W-Zin-We-''^ 



(10) 



In the limit of large /? no light is transmitted and all remains in the material during the period f < 0. 
Integrating Eq. [9] with respect to w gives an expression for Ea{z,t = 0) in the form of a convolution Fourier 
transforming along the spatial coordinate we get. 



Eaik,t = 0) = -/in_(--)sgn(fc)/3 
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Here fin{k/ri) is the input field at the time t — k/rj and r(^) is the gamma function. An expression for a at the time 
the field was fiipped can easily be obtained from the above result along with Eq. [21 

We now turn our attention to Region B and calculate ab and Ef, in that region subject to our desired boundary 
conditions ab{z,t = oo) = 0, Eb{z = zo,t) = /out+(^)- '^^^ output field /out+(^) ^^ ^^ ^^^^ stage undetermined. 

Solving Eq. [3] subject to the output condition ab{z,t = oo) = gives 
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{z,t) ^ ig dt'e''^''^'-*"^Eb{z,t') 
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dt' H{t' ~t)e-"i'^*''*')Eb{z,t') 



Fourier transforming gives 



abiz,Lu) = -igEb{z,Lu) 



{uj - -qz) 



+ ttS{uj — rjz) 



Substituting this in Eq. |4] and integrating gives 



(12) 



(13) 



Eb{z,uj) = Eb{z = ZQ,u!)exp dz' g'^N { + TrSiuj - -qz') 



-Pout+('^)exp(-/37r(iJ(cj - t]z) - H{uj - ■qzo))) 



u! — rjz 



Lo - r]Zo 
In a similar manner to in Region A, in the limit of large zq this can be approximated. 
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Eb{z,uj) = FQut^(cj)cxp(-/37riJ(cj - rjz)) 

Prom the above expression Eb{z — —zO,t) can be calculated resulting in 

fin+it) = /out+(i)e~'''^ 
Also from Eq.fTSlone can find and expression for Eb{k, t ~ 0): 
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Equations [TBI and [T7l give conditions, that if satisfied, will lead to a particular output pulse /out+(0 ^^"^ ^^^ the 
atoms being left in the ground state at the end of the process. From Eq. [I6]and comparing Eqns. [Tl]and[l7| we see 
that if an auxiliary input pulse is applied an output pulse related to our input pulse by 

foUt) = ~frA-m'^'^^y (18) 

Because |tp*^r(j/?)/r(— i/?) has a modulus of 1 for all t, it can be seen that the envelope of the output pulse is a 
time reversed version of the input. The required auxiliary pulse is given by [171 the practical usefulness of the memory 
would of course be greatly hampered by the need for this auxihary input pulse. Which must be an attenuated copy of 
the output pulse applied at the same time as the output pulse. However in the situation where the sample is optically 
thick, exp(— /37r) ^ 1, the required input pulse is zero. When arriving at Eq. [TSlone only needs to ensure that for 
t — the values for E match, the values for a will then also match because of Eqns. [2| and [4l 

The phase between the input pulse and the output pulse changes across the pulse an amount given by 
2/31og(tend/^start) where tstart and iend ai'G the start and end times times for the output pulse. This will be a modest 
phase shift in most situations where the memory might be used. A value for /? of 2 would provide sufficient optical 
depth for 99.999% efficiency, with this and 2 for iend/^start the phase shift across the pulse would be less than tt. This 
phase shift could be corrected for by some time dependent change in the optical path length, for instance a mirror 
mounted on a piezo or and electro optic modulator. Alternately two memories could be used in series with the initial 
frequency gradients opposite for each memory, the phase shifts of the two memories would then cancel. 

III. EFFICIENCY AND ANALOGUE BEAMSPLITTERS 

The fact that we can still get the desired output pulse with an optically thin sample by applying and auxiliary 
input pulse may not be all that helpful in the practical operation of the memory. It does however, along with the 
linearity of the differential equations, enable one to determine the effect of the memory in the optically thin regime. 
The difference between the true output pulse when working in the optically thin regime and the ideal case is equal 
to the effect of applying the auxiliary pulse alone. There is a close analogy between the operation of the memory 
in the finite optical depth regime and a beam splitter. An optical beam splitter takes a pairs of optical modes and 
outputs linear combinations of these modes. The action of the memory for both the times t < and i > can be 
reduced to that of a beam splitter, these are illustrated in Fig. [2] . Each of the two beamsplitters have two input and 
two output ports. For the left hand beamsplitter bottom input port, labeled /^"(i) and the mode of interest is the 
temporal mode of the optical wavepacket that we wish to store. The top, output port of the left hand beam splitter, 
labeled fouti^) for this represents the light that is transmitted through the sample, from Eq. [TOl we can see that this 
transmitted light has the same temporal mode as the input fight and that the amplitude transmitivity of our analogue 
beam splitter will be exp(— /^tt). The left hand input port and the right hand output ports of the two beamsplitters 
aren't optical modes but instead the beamsplitter acts on, and produces, 'polariton' excitations that have both an 
atomic and optical component. Like the polaritons considered in electromagnetically induced transparency (EIT) [17| 
these are combinations of photon and atom excitations. Unfike EIT polaritons optical gradient polaritons propagate 
in reciprocal space, to higher spatial frequencies, rather than propagating along the propagation direction [l4 |. 

The left hand port of the left hand beamsplitter represents the initial state of the optical field in the sample and the 
atoms before the input pulse is applied in for the case of our memory, during the operation of the memory this state 
win be a vacuum state. The right hand output mode of the first beamsplitter represents the mode of the excitation in 
the sample at the time the field is switched (t = 0) . The optical component of mode function is given by Eq. [TT] and 
the atomic part can be found from Eq. [2 One explanation for why the memory works is that output polariton mode 
for the left hand beamsplitter that represents the {t < 0) evolution matches the input polariton mode for righthand 
beamsplitter that represents the (t>0) evolution. So long as we choose input and output optical modes for the given 
by equations Eq. [16] and [HI The amplitude transmitivity of this second beamsplitter is given by exp(— /?7r). 

In order to have an efficient memory we would require excitation fed into the f^{t) port of the network shown 
in Fig. [2] to be directed entirely to the f^^^ (t) port . This requires high refiectivity for our analogue beamsplitters 
or equivalently large optical depth. The beamsplitter analogy easily enables calculation of the efficiency of the echo 
process as a function of optical depth. Because both beamsplitters have an amplitude transmitivity of exp(— /^tt) we 
end up efficiency for the echo given by 

EflSciency = (1 - exp{-2Pn)f 



This result agrees with the numerical simulations we have presented previously [13]. From the beamsplitter analogue 
one can also see immediately where the rest of the incident energy goes, exp(— 2/37r) was transmitted by the sample 
and exp(— 2/37r)(l — exp(— 2/37r)) remains in the sample. 
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Figure 2: The optical gradient echo memory represented as a pair of beam splitters. The bold arrows represent the flow of 
energy in the optically thick case. The spatial and temporal modes that correspond to the each port of the beam splitter is 
discussed in the text. 



IV. TRANSVERSE BROADENING AND MULTIPLE SWITCHING. 



Initial theoretical treatments of the efRciency of two-level controhed reversible inhomogeneous broadening (CRIB) 
echoes [lH investigated situations where the controlled detunings of the ions were not correlated with the position 
the atom. Such transverse broadening would arise from a microscopic broadening mechanism. This situation can be 
modeled by a large number of optically thin gradient echo sub- memories in series. The relevant network of analogue 
beamsplitters is shown in Fig. [31 If there are M sub memories, the transmitted input pulse will be attenuated by 
exp(— 2/37rM). There are M paths a photon can take from the input pulse f^{t) to the output pulse /Jut(^)) each of 
which involves two reflections and M — Itransmissions from our analogue beamsplitters. Because the echo from each 
sub-memory combines in phase with that from the previous memory these paths all combine constructively to give 
an efficiency for the echo of M^(l — exp(— 2/37r))^ cxp(— 2/37r(M — 1)). Taking the limit where each sub- memory is 
optically thin (/3 small) we arrive at an efficiency for the memory with microscopic broadening 

EfRciency = AP^n^M^ exp(-2/37rM) = d^g-d 

Here d = 2(3ttM is the optical depth, (the logarithm of the ratio of the energies of the incident an transmitted 
pulses). This result is the same as that derived in []J| by solution of the Maxwell-Bloch equations. 

Another situation to which our beamsplitter analogy can be applied is multiple switching of the broadening, both 
for longitudinal and transverse broadenings. If the broadening is switched in polarity after the first echo the excitation 
remaining in the sample will again be rephased leading to another echo. Multiple echoes of the original input pulse 
can be created in this way. The analogue beamsplitter network for multiple switching in a gradient echo is shown 
in Fig. m from this network one can see that the the fraction of the input energy coming out in the A;*'^echo is 
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Figure 3: Network of analogue beamsplitters relevant to a transverse CRIB memory modeled as a large number of optically 
gradient echo memories. 



(1 — exp(— 2/37r))^ exp(— 2/37r(fc — 1)). From the beamsplitter network it can be seen that with multiple swicthing the 
energy eventually leaves the sample in either as the transmitted pulse or as one of the echoes. The efRciencies of the 
echos as a function of /3 are plotted in Fig. [5l 

The beam splitter analogy can be extended further to the case of multiple switching from a system with microscopic 
broadening. The network of analogue beam splitters in this situation is similar to Figs. [3] but with a 2D array of 
beam splitters. The amplitude of each of the multiple echoes can be calculated by summing the amplitudes of all 
the possible paths through the beamsplitter network. This leads to the following expression for the portion of the 
incident energy that is output as the p^^ echo 
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Figure 4: Network of analogue beamsplitters relevant to an optical gradient echo memory where the field is switched a number 
of times 
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Here r = (1 — exp(— 2/37r) and i = cxp(— /?7r) are the amplitude transmission and reflection coefficients of the 
analogue beamsplitters. Taking the limit a large number of optically thin memories. We get the following expression 
for the efficiency of the p**" echo as 



Cp = exp(— d) 
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where d is the optical depth. As is shown in Fig. [6] combined efficiencies significantly greater than 54% can be achieved 
with multiple switching. 

V. CONCLUSION 

In conclusion we have presented an analytic treatment of the optical gradient echo, or longitudinal CRIB echo. We 
have shown that it is completely efficient in the case of large optical depth and that recall efficiencies of 100% can 
also be obtained for optically thin samples by multiple switching. By modehng a system with transverse CRIB by a 
large number of optically thin gradients can calculate echo efficiencies in this case also. As has been shown elsewhere 
the maximum echo efficiency is 4/e^ « 54%. Multiple switching can improve this overall efficiency with efficiencies of 
> 90% possible for the sum of the first 100 echoes. 
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Figure 5: (Color online) Echo efficiency for an optical gradient echo as a function of /3 showing the effect of multiple switching 
The ratio of the energies of the input and transmitted pulse is given by exp(— 27r/3). 
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Figure 6: (Color online) Echo efficiency as a function of optical depth for transverse broadening CRIB echo showing the effects of 
multiple switching. As is shown in the figure efficiencies significantly greater than 54% can be achieved with multiple switching 
but 100% efficiency is approached relatively slowly as the number of echoes increases. Approximately 100 echoes are needed to 
achieve greater than 90% efficiency. 



